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ABSTRACT 

We present the results from a Submillimeter Array survey of the 887 /um continuum 
emission from the protoplanetary disks around 95 young stars in the young cluster 
NGC 2024. Emission was detected from 22 infrared sources, with flux densities from 
~5 to 330 mJy; upper limits (at 3 a) for the other 73 sources range from 3 to 24 mJy. For 
standard assumptions, the corresponding disk masses range from ~0.003 to 0.2 Mq, with 
upper limits at 0.002-0.01 Mq. The NGC 2024 sample has a slightly more populated 
tail at the high end of its disk mass distribution compared to other clusters, but without 
more information on the nature of the sample hosts it remains unclear if this difference 
is statistically significant or a superficial selection effect. Unlike in the Orion Trapezium, 
there is no evidence for a disk mass dependence on the (projected) separation from the 
massive star IRS 2b in the NGC 2024 cluster. We suggest that this is due to either the 
cluster youth or a comparatively weaker photoionizing radiation field. 

Subject headings: circumstellar matter — planetary systems; protoplanetary disks — 
solar system: formation — stars: pre-main sequence 


1. Introduction 

The fundamental properties of circumstellar disks play a critical role in the formation and 
evolution of planets. While detailed knowledge of disk properties has come from extensive studies 
of nearby associations like Taurus-Auriga and p Ophiuchus (e.g., Beckwith et al. 1990; Osterloh & 
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Beckwith 1995; Andrews & Williams 2005, 2007; Andrews et al. 2009, 2010), stars in these regions 
form in loose agglomerations and are generally unaffected by their external environment. The 
majority of stars in the galaxy, including the Sun, formed in densely populated rich clusters (Lada 
&: Lada 2003; Porras et al. 2003; Williams 2010). The high stellar density and ultraviolet radiation 
from nearby massive stars in these regions can affect disk properties, threatening their development 
and potentially limiting their lifespans (Bonnell et al. 2003; Johnstone et al. 1998). Probing disk 
evolution in rich clusters is therefore crucial to our understanding of planet formation. 

The Orion star-forming complex contains the nearest rich clusters with massive stars, and is 
arguably the best region for studying how disk properties are affected by their environment. It is 
home to the clusters NGC 2024 (~0.5 Myr; Meyer 1996; Ali et al. 1998; Levine et al. 2006) and the 
Orion Nebula Cluster (ONC, ~l-2Myr; Reggiani et al. 2011; Da Rio et al. 2010). Near-infrared 
observations of the NGC 2024 cluster members have revealed 233 young stars (Meyer 1996), of which 
> 85% exhibit an infrared excess indicative of warm dust in the inner regions of protoplanetary 
disks (Haisch et al. 2000, 2001). These stars are still deeply embedded in molecular cloud material 
(Barnes et al. 1989; Lada 1991), in line with their suggested extreme youth. The most massive star 
in the region is thought to be IRS 2b, with a spectral type in the range of 08 to B2 (Bik et al. 
2003; Barnes et al. 1989). Although the earlier end of that range could more easily explain the 
radio continuum emission in the region, it has been suggested that a collection of slightly lower 
mass stars could together be responsible for the total ionizing flux (Meyer et al. 2008). IRS 2b is 
located 5" northwest of the early B-type star IRS 2 (Grasdalen 1974), the brightest infrared and 
radio source in NGC 2024 (Barnes et al. 1989; Rodriguez et al. 2003). 

The relative youth of NGC 2024 makes it a particularly appealing region to probe the initial 
properties of disks. Through comparisons with the older ONC (and similar clusters, e.g. a Ori; 
Williams et al. 2013), we can constrain key timescales for disk evolution in rich clusters. With 
reference to pre-main sequence evolution models, the positions of the NGC 2024 members in a color- 
magnitude diagram indicate very young ages: the Baraffe et al. (1998) models suggest <1 Myr, and 
the D’Antona Sz Mazzitelli (1997) models argue for ~0.5Myr (Meyer 1996; Ali et al. 1998; Levine 
et al. 2006). The absolute ages of young stars are highly uncertain (e.g., Soderblom et al. 2013). 
However, in a relative sense these are significantly younger ages than have been estimated for other 
clusters using the same technique and models (e.g., see Eisner & Carpenter 2003). They overlap 
with the earliest stages of the evolution process for circumstellar material (e.g., Evans et al. 2009). 

Although the infrared excess emission found for most NGC 2024 members confirms the ubiquity 
of disks in this region, it traces only a small fraction of the disk material. The low optical depths for 
the continuum emission at longer, (sub)millimeter wavelengths is required to quantitatively probe 
the masses of these disks. While the disk population in the neighboring ONC has now been studied 
in some detail at (sub)millimeter wavelengths with interferometers (Williams et al. 2005; Eisner & 
Carpenter 2006; Eisner et al. 2008; Mann & Williams 2010; Mann et al. 2014), the disks in NGC 
2024 are relatively unexplored. Eisner &: Carpenter (2003) made the sole attempt at measuring the 
NGC 2024 disk mass distribution, using a 3 mm survey of 150 targets with modest sensitivity (a disk 



Fig. 1.— (left) A 10'x 10' false-color near-infrared image of NGC 2024, also known as the Flame 
Nebula (Meyer et al. 2008). The white box shows the region targeted in the SMA survey, (right) 
A JCMT-SCUBA 850 /rm image of the inset region (obtained from Di Francesco et al. 2008), with 
the pointing locations and dimensions of the SMA primary beam overlaid as white circles, labeled 
as in Table 1. Crosses show the location of young stars identified in A-band imaging (Meyer 1996). 
The most massive star of the cluster, IRS 2b, is labeled in both panels. 


mass upper limit of ~0.035Mq). They detected two massive disks (0.08 and 0.24 Mq; signihcantly 
larger than seen in the ONC), and argued that image-stacking suggested that the average disk mass 
was ~0.005 Mq (comparable to Taurus-Auriga and p Ophiuchus). 

Here we present the results of a new Submillimeter Array (SMA) survey of the 887 ^m con¬ 
tinuum emission toward 95 young stars in the NGC 2024 cluster. Leveraging the steep scaling 
between the continuum flux and observing frequency (Fy oc these observations represent an 

order of magnitude improvement in sensitivity over the previous work in this region, and should be 
capable of detecting the average disk mass as suggested by Eisner &: Carpenter (2003). This survey 
represents the deepest attempt to measure the disk mass distribution in a very young, rich cluster, 
and thereby to probe how environment impacts basic disk properties. The observations and their 
calibration are described in Section 2. The flux measurements and their estimated conversion to 
disk masses are presented in Section 3. We make a comparison of the derived disk mass distribution 
with other regions, examine the dependence of disk mass on location in the cluster, and discuss the 
implications for planet formation in rich clusters in Section 4. 


2. Observations 

Observations of nine distinct pointings containing a total of 95 young stars were conducted 
with the SMA (Ho et al. 2004) in the fall of 2011, using the compact array configuration (baselines 
of ~8-50m). The pointing centers are listed in Table 1 and shown in Figure 1, and were chosen 
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Fig. 2.— Synthesized images of the 887 /xm continuum emission toward the 22 detected infrared 
sources (see Table 2). Each panel is 12" (~5000AU) on a side, and includes labels in the Meyer 
(1996) IRC designation. Contours are drawn at 3 a intervals, and the synthesized beam dimensions 
are shown in the bottom right corners of each panel. In all but three panels, the image center 
corresponds to the infrared source position; otherwise, those positions are marked with crosses. 
The third panel includes the two cases with multiple possible identifications; see Table 2. 

to maximize the number of young stars imaged in each (35" FWHM) primary beam while mini¬ 
mizing contamination from the bright, non-uniform molecular cloud background in the region (see 
Section 3). The SMA double sideband receivers were tuned to a local oscillator (LO) frequency 
of 338.213 GHz (887/xm, see Table 1). Each sideband provided 4 GHz of bandwidth, centered ±5 
GHz from the LO frequency. Three observing tracks were shared between three separate pointings 
in each. The observations of the NGC 2024 fields were interleaved with nearby gain calibrators on 
15 minute intervals. Weather conditions for all observations were good, with <2 mm of precipitable 
water vapor, r(225 GHz) <0.1. Table 1 summarizes the relevant observational information. 

The raw visibilities were calibrated using the MIR software package. Passband calibration was 
conducted using the bright, compact radio sources, 3G 279, 3G 84, or J0854-I-201. The absolute 
flux scale was derived from observations of Titan and Uranus (and checked against 3G 111), and 
is accurate to ~10%. Amplitude and phase calibration were performed using observations of the 
nearby sources J0423—013, J0530-I-135, and J0607—085. The calibrated visibilities were naturally 
weighted and Fourier inverted, then CLEANed to generate the synthesized continuum maps shown 
in Figure 2 using MIRIAD (Sault et al. 1995). The synthesized maps were created after eliminating 
projected antenna spacings shorter than 27 kA, to filter out extended emission on size scales > 7.5". 













- 5 - 


This scale was chosen to preserve the compact emission from the disks, while minimizing contami¬ 
nation from the bright, extended cloud background; it is the same as used in the analysis of dust 
emission from disks in the ONC (Mann &: Williams 2010) for the sake of consistency. Simulations 
of the background (see Section 3) confirmed that the 27 kA cutoff resolves out most of the extended 
emission and thereby helps reduce the effective RMS noise levels by roughly a factor of ~2. 


3. Results 

Continuum emission was detected toward 22 near-infrared sources at >3x the measured RMS 
noise level (see Table 2).^ The observed flux density. Fobs, for each source was determined with a 
Gaussian fit in the image plane, after correcting for primary beam attenuation. The detected sources 
were associated with their stellar counterparts by reference to the near-infrared catalog of Meyer 
(1996). In two cases, IRC 106/098 and IRC 101/093 (see Table 2), it is difficult to unambiguously 
distinguish two potential associations with infrared sources. For the former, observations with 
improved sensitivity and resolution (and ideally Fourier sampling) will be required to robustly 
differentiate the options. On the other hand, it is entirely possible that the very bright emission 
near IRC 101/093 is unassociated with either source, but instead is tracing an embedded source 
at an earlier evolutionary stage. Another 73 infrared sources were covered in the survey, but 
not detected. Limits on their sub millimeter continuum flux densities were determined from local 
measurements of the RMS noise level, and are also listed in Table 2. 

The observed flux densities are the linear combination of several emission contributors; free- 
free radiation from ionized material (Fg), and thermal radiation from dust in the surrounding 
molecular cloud (F/ioud) and the disk (Tdisk)- We assumed that any free-free emission is optically 
thin, with a spectrum Fg oc and used deep VLA 3.6 cm measurements (Rodriguez et al. 

2003, see Table 2) when available as normalizations to extrapolate Fg up to 887 fiui. Only four 
of the VLA 3.6 cm sources overlap with the SMA 887 fiui detections. Upper limits at 3.6 cm are 
sufficiently low (~50 /rJy at 3 a) that extrapolated estimates of free-free contamination for the other 
targets are considered negligible. The observed targets are embedded in their host cloud, which 
itself produces significant dust emission on large spatial scales. To estimate that cloud contribution 
at each location in the SMA maps, we simulate the SMA response to its large-scale emission as 
observed with the SCUBA instrument on the James Clerk Maxwell Telescope (JCMT; Di Francesco 
et al. 2008), following the approach of Mann &: Williams (2010). For each SMA pointing, we Fourier 
transformed the appropriate JCMT map, sampled the emission onto the observed u, u-tracks, and 
generated CLEANed maps. The cloud emission, F/ioud, was then estimated toward each source 


^An additional ~20 emission peaks (each at ~3 (j) were identified, but are not coincident with any known near- 
infrared sources (some of them are visible in Fig. 2). This is the expected number of 3 cr noise peaks over such a large 
survey area (perhaps even an under-estimate when considering the sparse Fourier sampling), although we cannot rule 
out the possibility that some are real sources associated with deeply embedded objects. 



location. We found that cloud contamination was typically <10%, considerably lower than in 
the ONC. After accounting for any free-free or cloud contamination, the remainder of F^hs was 
associated with the disk, Fdisk- Table 2 lists the decomposed contributions for each source. 


Relatively little is known about the nature of the infrared sources associated with submillimeter 
emission in this survey. Given their neutral or modestly red JHK colors, we are making the 
assumption that the vast majority of them are in the so-called Class II or T Tauri stage, young 
stars with disks but no remnant envelope material. It is possible that the few most luminous 
sources could be in the earlier Class I stage, if their observed near-infrared emission is primarily 
tracing scattered light from outflow cavities in their envelope structures (Eisner 2012; Sheehan &: 
Eisner 2014). As we noted above, the very brightest source in this sample cannot be unambiguously 
associated with an infrared source; it is plausible that the emission may originate in a dense envelope 
around a Class 0 protostar. Without a more complete set of ancillary information (e.g., full spectral 
energy distributions), a refined classification is not yet feasible. 


Eollowing Eisner &: Carpenter (2003), we aimed to constrain the mean emission level of the 73 
undetected sources in the SMA fields by “stacking” the data. However, there is sufficient concern 
with combining data in the image plane when the individual fields were sampled so sparsely in 
the Eourier domain. Instead we performed a complementary analysis on the visibilities. First, we 
removed the detected sources by subtracting Gaussian models of the emission from the observed 
visibilities. We then generated 73 permutations of the visibility data, each with a phase shift that 
accounts for the location of the undetected cluster member.^ Those permutations were co-added 
and then imaged as described in Section 2. We found no emission associated with this stacked 
dataset, and placed a 3cr upper limit of ~2.4mJy on the ensemble average. Assuming a typical 
spectrum that scales like i% oc v‘^~^, this limit is ~4~I0x lower than the ensemble mean flux density 
that was estimated from image-plane stacking by Eisner &: Carpenter (2003). 


We assume that the observed continuum emission is (mostly) optically thin, and therefore a 
sensitive probe of the dust mass. Since most of the emission originates in the cool, nearly isothermal 
outer regions of a disk, we estimate the mass as 


Tldisk 


d Edisk 


( 1 ) 


where d is the distance, is the opacity per gram of the disk material, and By{T) is the Planck 
function at a characteristic temperature (e.g., Beckwith et al. 1990). Disk masses were calculated for 
the 22 detected sources (see Table 2) using Eq. I and standard assumptions (for ease of comparison 
with other studies): a characteristic dust temperature T = 20K, the Beckwith et al. (1990) opacity 
Kjy = 0.034cm^g“^ at 887/rm (which implicitly assumes a 100:1 gas-to-dust mass ratio), and a 
distance d = 415 pc to NGC 2024, based on observations of B-type stars in the cluster (Anthony- 
Twarog 1982). Disk emission was inferred toward the source IRS 2 (IRC 232), an early B-type star; 


^Technically this includes substantial duplication of the data. However, the individual sources are located far 
enough apart that their mutual contributions at any given phase shift are considered negligible. 
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in calculating its mass, we adopted a higher dust temperature of 40 K (e.g., see Beuther et al. 2002; 
Sridharan et al. 2002). Systematic uncertainties in the Mdisk estimates are dominated by the poorly 
constrained values of the dust opacities, which are ambiguous at the order of magnitude level (e.g., 
Henning Sz Stognienko 1996). Some relatively minor additional uncertainties from optical depth 
effects could also contribute, especially given the unknown disk structures: Andrews Sz Williams 
(2005) estimated that these are on the order of ~10% for the fainter luminosities that characterize 
this sample, but could rise to as high as ~50% at the high luminosity end. 

The overall disk mass sensitivity of this SMA survey depends on the locations of each target 
relative to the field (pointing) center, the varying levels of cloud emission, and any free-free emission 
contributions. A mass completeness level for the survey was estimated based on Monte Carlo 
simulations. Synthetic disks (point sources) with emission appropriate for a given Mdisk (see Eq. 1) 
were injected into the large-scale emission maps from the JCMT, Fourier inverted and sampled onto 
the observed spatial frequencies, and then each field was imaged as in Section 2. By measuring the 
fraction of synthetic targets that were detected (>3 (t) for each input Mdisk in these Monte Carlo 
simulations, we determined that the survey is essentially 100% complete for Mdisk > 0.01 Mq, 
and roughly 50% complete for Mdisk > 0.004 Mq. The upper limit on the stacked ensemble of 
undetected sources described above corresponds to ~0.001 Mq. 

Four of the 13 sources detected by Eisner & Carpenter (2003) at 3 mm were also detected in 
the SMA survey at 887/rm: IRC 124, IRC 071, IRC 101, and IRS 2 (their sources 1, 4, 8, and 9, 
respectively). The remaining nine sources were either outside our survey area (their sources 3, 11, 
12, and 13) or were not detected with the SMA (IDs 2, 5, 6, 7, 9, and 10); the latter cases likely 
indicate that the radiation detected at 3 mm was generated by a non-dust emission mechanism. The 
observed flux ratios can be used to constrain the spectral slopes between 887 /um and 3 mm of the 
targets with overlapping detections. We hnd spectral indices, a where i% oc zz“, of 2.8T0.2 (IRC 124 
= source 1), 2.8 ± 0.1 (IRC 071 = source 4), 2.6 ± 0.1 (IRC 101 = source 8), and 0.1 ± 0.1 (IRS 2 
= source 9); the IRS 2 measurement is obviously strongly impacted by free-free contamination at 
3 mm. These spectral indices are at the higher end of the distributions inferred in ~l-3Myr-old 
star-forming clusters (Ricci et al. 2010b,a, 2011a, b), perhaps hinting at some age evolution in the 
disk-integrated opacity spectrum due to dust grain growth (e.g., Miyake &: Nakagawa 1993; Draine 
2006). More measurements for young sources in NGC 2024 would be desireable. 


4. Discussion 

We conducted a large-scale survey of the 887 ^m continuum emission toward 95 young stars in 
the ~0.3Myr-old NGC 2024 cluster using the SMA. Assuming standard conversions for optically 
thin dust emission, this survey is complete down to a disk mass limit of ~0.01 Mq (3 a), although 
in some regions the survey is slightly more sensitive. We detected a total of 22 disks (23 ± 5% of 
the sample), including four that were previously detected at 3mm by Eisner &: Carpenter (2003). 
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Fig. 3.— The cumulative distribution of disk masses in the NGC 2024 survey, constructed with 
the Kaplan-Meier product limit estimator to account for upper limits. The estimated 50 and 100% 
completeness levels are marked as dotted vertical lines. The analogous distribution for the Taurus 
star-forming region (Andrews et al. 2013) is shown in gray for reference. The upper end of the 
NGC 2024 disk mass distribution favors slightly higher masses, as might be expected for a younger 
cluster, although the selection effects of this sample are not yet well characterized. 

Figure 3 shows the cumulative distribution of disk masses in this survey, where we have in¬ 
corporated the upper limits by employing the Kaplan-Meier product limit estimator for a censored 
sample (Feigelson & Nelson 1985). Although these data are sensitivity-limited to probe only rel¬ 
atively massive disks, we find that the fraction of disks with large Mdisk is relatively high: ~20% 
have Mdisk > 0.01 M©, and ~10% have Mdisk > 0.1 Mq. Taken at face value, this suggests that 
the high-mass tail of the NGC 2024 Mdisk distribution is more populated than in slightly older 
(~l-3Myr) clusters like the ONC (Mann et al. 2014), Ophiuchus (Andrews & Williams 2007), and 
Taurus (see the corresponding distribution function in Fig. 3 for a direct comparison; Andrews & 
Williams 2005; Andrews et al. 2013). The standard censored two-sample tests advocated by Feigel¬ 
son &: Nelson (1985) indicate a marginal (~2 (t) quantitative offset, with the NGC 2024 distribution 
shifted to higher masses by a factor of ~ 1.5-2. 

However, such comparisons can be misleading if they do not account for selection biases. Disk 
masses are known to depend on factors like the (host) mass (Andrews et al. 2013), multiplicity 
(Harris et al. 2012; Akeson &: Jensen 2014), and evolutionary state (Andrews & Williams 2005) of 
the target. Unfortunately, little is known about these properties for the NGC 2024 sample. If we 
make the assumptions that this survey has targets drawn from the same host mass function and 
with the same multiplicity statistics as in Taurus, and suggest that all of the near-infrared sources 
in the SMA fields harbor disks with no envelopes, we can use the Monte Carlo approach advocated 
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by Andrews et al. (2013) to compare with the reference M^isk distribution in Taurus. The results 
indicate that the two M^isk distributions are statistically indistinguishable (given the current data), 
suggesting that there is relatively little evolution in Mdisk up to a few Myr. More robust constraints 
on changes in this distribution would require a better characterization of the NGC 2024 targets, 
and secondarily an expanded continuum census. 

But independent of these (potential) selection effects, a particularly interesting comparison can 
still be made between the disks in the NGC 2024 cluster and the ONC, since both regions host high- 
mass stars in their immediate environments that could potentially modify the disk mass distribution. 
Sources in the ONC found within ~0.03pc of the massive star 9^ Ori C have systematically lower 
Mdisk than those at larger separations (Mann &: Williams 2009, 2010; Mann et al. 2014), reflecting 
the consequences of external evaporation from strong photoionizing sources on disk dissipation 
timescales (Johnstone et al. 1998; Storzer & Hollenbach 1999; Richling & Yorke 2000; Scally &: 
Clarke 2001; Matsuyama et al. 2003; Adams et al. 2004). Figure 4 shows Mjisk as a function of 
the projected distance from the most massive stars in both the NGC 2024 cluster and the ONC, 
IRS 2b and 0^ Ori C, respectively. Unlike the ONC, we find no evidence of a distance-dependent 
disk mass distribution in NGC 2024. Although the total number of disk detections in the NGC 
2024 region is limited, several massive disks identified here are located <0.01 pc from IRS 2b. 
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Fig. 4.— Disk masses in the NGC 2024 (black) and ONC (gray; see Mann &: Williams 2010; 
Mann et al. 2014) clusters as a function of their projected separations from the nearest massive 
star, IRS 2b and 6^ Ori C, respectively. Circles represent submillimeter continuum detections of 
dust disk emission, and horizontal line segments mark 3 a upper limits. The depletion in M^isk at 
small projected separations seen for the ONC is not apparent for the NGC 2024 cluster, perhaps 
because it represents an earlier evolutionary stage or due to the (presumably) comparatively weaker 
photoionizing radiation field present. 
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We consider two likely, and not mutually exclusive, reasons for this difference between NGC 
2024 and the ONC. First is an evolutionary argument, based on the fact that NGC 2024 appears to 
be considerably younger than the ONC (see Meyer 1996; Eisner &: Carpenter 2003). That relative 
youth could mean that we are observing more of a primordial M^isk distribution in NGC 2024, 
before external evolutionary processes like photoevaporation have had time to make a significant 
impact. Second is an environmental distinction, in that IRS 2b is an intrinsically less luminous 
source than 6^ Ori C, and therefore could produce a substantially weaker photoionizing radiation 
field that is less capable of stripping material from its surrounding disks. Bik et al. (2003) suggested 
that IRS 2b has a spectral type of 08-B2, although the radio continuum flux measured by Barnes 
et al. (1989) indicates that the earlier type is more appropriate. However, Meyer et al. (1997) 
instead suggested that early B-types are the preferred spectroscopic classification for IRS 2b, and 
that it may not be the sole or even dominant source of ionizing radiation in the region (see also 
Meyer et al. 2008). A small group of weaker ionizing sources might provide sufficiently attenuated 
mass-loss rates and explain the absence of a separation dependence on the Mdisk distribution. Even 
if IRS 2b was more luminous than described, NGC 2024 is more heavily extincted than the ONC, 
and the high energy photons could be easily absorbed by the cloud. 

Euture, more sensitive observations of NGC 2024 with the Atacama Large Millimeter Array 
(ALMA) will permit an exploration of the full disk mass distribution in this young region, facili¬ 
tating stronger constraints on disk dissipation by ultraviolet photoevaporation, on the evolutionary 
timescales of disks in rich clusters, and on the initial conditions of the planet formation process 
for the majority of stars in the galaxy. Such campaigns should necessarily be coupled with a more 
comprehensive characterization of the NGC 2024 stellar population. 


We thank the referee for a very helpful review. The Submillimeter Array is a joint project 
between the Submillimeter Astrophysical Observatory and the Academica Sinica Institute of As¬ 
tronomy and Astrophysics and is funded by the Smithsonian Institution and the Academica Sinica. 
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Table 1. Summary of Submillimeter Array Observations 


Field 

a (J2000) 

(J2000) 

UT Date 

T 

a (mJy/beam) 

Ob (") 

PAfe (°) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

A 

05 41 47.0 

-01 54 47 

2011 Sep 23 

0.06-0.07 

1.0 

2.5 X 1.5 

158 

B 

05 41 45.4 

-01 54 18 

2011 Sep 23 

0.06-0.07 

3.0 

2.5 X 1.5 

158 

C 

05 41 45.1 

-01 54 54 

2011 Sep 23 

0.06-0.07 

3.5 

2.5 X 1.5 

158 

D 

05 41 45.3 

-01 53 47 

2011 Oct 18 

0.03-0.07 

1.4 

1.8 X 1.5 

15 

E 

05 41 43.5 

-01 53 30 

2011 Oct 18 

0.03-0.07 

1.0 

1.8 X 1.5 

15 

F 

05 41 43.2 

-01 54 04 

2011 Oct 18 

0.03-0.07 

4.6 

1.8 X 1.5 

15 

G 

05 41 41.2 

-01 53 40 

2011 Dec 29 

0.08-0.15 

2.2 

2.2 X 1.7 

5 

H 

05 41 39.2 

-01 54 02 

2011 Dec 29 

0.08-0.15 

1.0 

2.2 X 1.7 

5 

I 

05 41 38.7 

-01 53 23 

2011 Dec 29 

0.08-0.15 

1.4 

2.2 X 1.7 

5 


Note. — All observations were conducted at 887 /rm using the compact array configuration. 
Col. (1): SMA pointing, as labeled in Figure 1. Cols. (2, 3): Phase center coordinates. Col. (4): 
UT Date of observation. Col. (6): Range of zenith optical depths at 225GHz. Col. (8); RMS 
noise level measured in emission-free regions within the primary beam of the naturally-weighted 
synthesized maps. Cols. (9, 10): Dimensions and orientations of the synthesized beams. 
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Table 2. Inferred Disk Fluxes and Masses 


Source 

Field 

a (J2000) 

(5 (J2000) 

.fobs 

-fff 

f^cloud 

f/iisk 

-^disk 

(IRC) 




(mJy) 

(mJy) 

(mJy) 

(mJy) 

(0.01 Mq) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

065 

A 

05 41 46.9 

-01 54 46.9 

7.2 ±1.2 

0.8 

0.5 

5.9 

0.3±0.1 

128 

B 

05 41 45.1 

-01 54 06.9 

21.2 ±3.0 


0 

21.8 

1.2±0.2 

106/098 

B 

05 41 45.6 

-01 54 22.4 

10.7 ±3.0 


0.3 

10.5 

0.6±0.2 

101/093* 

B 

05 41 45.4 

-01 54 26.3 

298.6 ±3.0 

0.7 

0.3 

297.7 

16.0±0.2 

232* (IRS 2) 

B 

05 41 45.8 

-01 54 30.2 

209.4 ± 3.0 

12.0 

0 

199.3 

4.0±0.2 

050 

C 

05 41 45.8 

-01 55 06.2 

50.7 ±3.0 


0.1 

50.6 

2.7±0.3 

067 

C 

05 41 45.1 

-01 54 47.3 

119.0 ±3.0 


0.2 

118.9 

6.4±0.2 

071* 

C 

05 41 44.1 

-01 54 45.8 

330.5 ±3.0 

0.2 

0.1 

330.1 

17.7±0.3 

229 

C 

05 41 44.0 

-01 54 43.1 

192.6 ±3.0 


0.2 

192.5 

10.3±0.4 

058 

C 

05 41 44.1 

-01 54 54.6 

34.1 ±3.0 


0.1 

34.0 

1.8±0.3 

047 

C 

05 41 44.1 

-01 55 06.5 

210.7 ±3.0 


0.1 

210.6 

11.3±0.4 

153 

D 

05 41 44.7 

-01 53 48.7 

55.0 ±1.2 


0 

55.2 

3.0±0.1 

134 

D 

05 41 44.7 

-01 54 01.6 

17.1 ± 1.2 


0.1 

17.0 

0.9±0.1 

197 

E 

05 41 43.5 

-01 53 24.8 

6.1 ± 1.0 


0.8 

5.3 

0.3±0.1 

168 

E 

05 41 44.1 

-01 53 42.3 

9.7 ± 1.0 


0.9 

8.8 

0.5±0.1 

117 

F 

05 41 43.2 

-01 54 15.6 

58.1 ±4.5 


0.5 

57.7 

3.1±0.3 

124* 

F 

05 41 42.5 

-01 54 08.8 

305.5 ±4.5 


0.9 

304.6 

16.4±0.3 

151 

G 

05 41 41.5 

-01 53 48.6 

179.5 ±2.0 


0.1 

179.4 

9.6±0.1 

179 

G 

05 41 41.7 

-01 53 35.2 

46.3 ±2.0 


0.0 

46.4 

2.5±0.1 

191 

G 

05 41 41.5 

-01 53 26.7 

32.3 ±2.0 


0 

33.1 

1.8±0.2 

182 

G 

05 41 41.1 

-01 53 31.2 

172.3 ±2.0 


0 

172.7 

9.3±0.1 

129 

H 

05 41 40.2 

-01 54 06.6 

6.3 ± 1.3 


0 

6.3 

0.3±0.1 

Non Detections 




3cj Upp Limit 





066 

A 

05 41 46.5 

-01 54 47.2 

4.0 

0.8 

0.1 

3.1 

0.2 

068 

A 

05 41 46.7 

-01 54 45.1 

3.8 


0.2 

3.6 

0.2 

062 

A 

05 41 46.8 

-01 54 49.3 

3.8 


0.2 

3.6 

0.2 

060 

A 

05 41 47.0 

-01 54 53.1 

3.9 


0.3 

3.6 

0.2 

064 

A 

05 41 47.3 

-01 54 47.1 

3.7 


0.1 

3.6 

0.2 

107 

B 

05 41 44.7 

-01 54 20.7 

12.0 


0.2 

11.8 

0.6 

231 

B 

05 41 44.7 

-01 54 31.4 

17.1 


0.2 

16.9 

0.9 

099 

B 

05 41 44.8 

-01 54 25.4 

12.0 

0.2 

0.1 

11.8 

0.6 
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Table 2—Continued 


Source 

Field 

a (J2000) 

5 (J2000) 

Fobs 

Fff 

Fcloud 

Fdisk 

F7(jisk 

(IRC) 




(mJy) 

(mJy) 

(mJy) 

(mJy) 

( 0.01 Mq ) 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

134 

B 

05 41 44.9 

-01 54 01.8 

17.8 


0.1 

17.7 

1.0 

102 

B 

05 41 45.2 

-01 54 23.1 

9.7 


0.2 

9.5 

0.5 

127 

B 

05 41 45.7 

-01 54 06.5 

12.7 


0.1 

12.6 

0.7 

122 

B 

05 41 45.9 

-01 54 11.3 

11.0 

0.1 

0.2 

10.7 

0.6 

114 

B 

05 41 46.1 

-01 54 15.0 

11.5 


0.1 

11.4 

0.6 

109 

B 

05 41 46.4 

-01 54 18.4 

14.2 


0.1 

14.0 

0.8 

073 

C 

05 41 44.9 

-01 54 42.6 

12.0 


0.1 

11.9 

0.6 

059 

C 

05 41 45.0 

-01 54 55.1 

9.1 

0.2 

0.2 

8.7 

0.5 

052 

C 

05 41 45.3 

-01 55 03.0 

10.9 


0.1 

10.8 

0.6 

063 

C 

05 41 45.7 

-01 54 49.7 

10.9 


0.2 

10.7 

0.6 

057 

C 

05 41 45.7 

-01 54 58.3 

11.0 

0.2 

0.1 

10.6 

0.6 

053 

C 

05 41 45.9 

-01 55 02.2 

11.3 


0.2 

14.2 

0.8 

157 

D 

05 41 44.1 

-01 53 47.5 

7.1 


0 

7.1 

0.4 

175 

D 

05 41 44.4 

-01 53 36.2 

7.0 


0.1 

6.9 

0.4 

138 

D 

05 41 44.5 

-01 53 57.4 

6.2 


0 

6.2 

0.3 

176 

D 

05 41 45.3 

-01 53 34.8 

4.9 


0.2 

4.7 

0.3 

139 

D 

05 41 45.4 

-01 53 56.5 

4.4 


0 

4.4 

0.2 

161 

D 

05 41 45.9 

-01 53 45.1 

4.2 


0.1 

4.1 

0.2 

167 

D 

05 41 46.1 

-01 53 42.3 

5.0 


0.1 

4.9 

0.3 

156 

D 

05 41 46.2 

-01 53 47.0 

5.4 


0.1 

5.3 

0.3 

206 

E 

05 41 42.9 

-01 53 16.3 

5.4 


0.1 

5.3 

0.3 

202 

E 

05 41 43.2 

-01 53 18.4 

4.3 


0.2 

4.1 

0.2 

181 

E 

05 41 43.3 

-01 53 32.8 

3.1 


0.2 

2.9 

0.2 

187 

E 

05 41 43.3 

-01 53 30.2 

3.0 


0.2 

2.8 

0.2 

193 

E 

05 41 43.5 

-01 53 26.4 

3.1 


0.1 

3.0 

0.2 

170 

E 

05 41 43.8 

-01 53 38.6 

3.7 


0.1 

3.6 

0.2 

188 

E 

05 41 43.9 

-01 53 29.3 

3.3 


0.1 

3.2 

0.2 

172 

E 

05 41 44.1 

-01 53 37.6 

4.1 


0 

4.1 

0.2 

175 

E 

05 41 44.4 

-01 53 36.2 

4.7 


0.1 

4.6 

0.2 

133 

E 

05 41 42.5 

-01 54 03.1 

17.1 


0.2 

16.9 

0.9 

115 

E 

05 41 42.8 

-01 54 16.1 

20.4 


0.1 

20.3 

1.0 
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Table 2—Continued 


Source 

Field 

a (J2000) 

5 (J2000) 

Fobs 

Fff 

Fcloud 

Fdisk 

FTdisk 

(IRC) 




(mJy) 

(mJy) 

(mJy) 

(mJy) 

(0.01 Mq) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

140 

F 

05 41 43.1 

-01 53 57.5 

14.9 


0.3 

14.6 

0.8 

121 

F 

05 41 43.4 

-01 54 12.5 

16.1 


0 

16.1 

0.9 

143 

F 

05 41 43.5 

-01 53 56.9 

15.9 

0.1 

0.2 

15.6 

0.8 

150 

F 

05 41 43.9 

-01 53 51.5 

23.9 


0.2 

23.7 

1.3 

132 

F 

05 41 44.0 

-01 54 03.2 

18.3 


0.3 

18.0 

1.0 

131 

F 

05 41 44.2 

-01 54 21.9 

21.9 


0.3 

21.6 

1.2 

178 

G 

05 41 40.0 

-01 53 35.7 

11.9 


0.1 

11.8 

0.6 

180 

G 

05 41 40.2 

-01 53 34.3 

10.8 


0.1 

10.7 

0.6 

194 

G 

05 41 40.5 

-01 53 27.4 

11.0 


0 

11.0 

0.6 

158 

G 

05 41 40.6 

-01 53 48.5 

8.3 


0 

8.3 

0.4 

165 

G 

05 41 41.7 

-01 53 44.6 

7.0 


0 

7.0 

0.4 

166 

G 

05 41 42.0 

-01 53 44.3 

8.4 


0.1 

8.3 

0.4 

185 

G 

05 41 42.1 

-01 53 32.1 

10.0 


0.1 

9.9 

0.5 

146 

H 

05 41 38.0 

-01 53 57.6 

7.8 


0.1 

7.7 

0.4 

148 

H 

05 41 38.9 

-01 53 54.0 

4.7 


0 

4.7 

0.2 

126 

H 

05 41 39.0 

-01 54 09.7 

4.5 


0 

4.5 

0.2 

145 

H 

05 41 39.1 

-01 53 58.6 

4.0 


0.1 

3.9 

0.2 

120 

H 

05 41 39.1 

-01 54 14.4 

5.4 


0.1 

5.3 

0.3 

136 

H 

05 41 39.2 

-01 54 02.4 

3.9 


0 

3.9 

0.2 

147 

H 

05 41 40.1 

-01 53 56.1 

6.3 


0.1 

6.2 

0.3 

213 

I 

05 41 37.9 

-01 53 11.5 

7.2 


0 

7.2 

0.4 

203 

I 

05 41 38.0 

-01 53 19.0 

5.3 


0.1 

5.2 

0.3 

198 

I 

05 41 38.1 

-01 53 25.8 

4.7 


0.1 

4.6 

0.2 

192 

I 

05 41 38.3 

-01 53 28.7 

4.5 


0.1 

4.4 

0.2 

186 

I 

05 41 38.3 

-01 53 33.5 

5.4 


0 

5.4 

0.3 

204 

I 

05 41 38.4 

-01 53 18.2 

4.3 


0 

4.3 

0.2 

200 

I 

05 41 38.6 

-01 53 23.1 

3.9 


0 

3.9 

0.2 

201 

I 

05 41 39.3 

-01 53 22.3 

4.7 


0.1 

4.6 

0.2 

195 

I 

05 41 39.5 

-01 53 27.1 

5.3 


0.1 

5.2 

0.3 

183 

I 

05 41 39.5 

-01 53 34.0 

6.7 


0 

6.7 

0.4 
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Note. — Sources associated with asterisks were also detected at 3 mm by Eisner h Car¬ 
penter (2003): IRC 106/098 = Source 8, IRC 232 = Source 9 (IRS2), IRC071 = Source 
4, and IRC 124 = Source 1. Col. (1): Source designation, according to Meyer (1996). 
Col. (2): SMA field, as labeled in Figure 1 and Table 1. Cols. (3, 4): SMA emission cen¬ 
troid coordinates. Col. (5): Integrated continuum flux density, or 3 cj upper limits on the 
non-detections, corrected for SMA primary beam attenuation. Col. (6): Extrapolated con¬ 
tribution of free-free emission at 887 /rm estimated from Rodriguez et al. (2003) measure¬ 
ments. Col. (7): Estimated contribution from large-scale cloud emission. Col. (8): Derived 
dust continuum flux density from the disk. Col. (9): Inferred disk mass (uncertainty does 
not include systematics in the absolute flux scale, which contribute an additional ~10%). 
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